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Abstract 

We discuss a grand unified model based on SUSY SU(5) in extra dimensions and 
on the flavour group ^4 x U(l) which, besides reproducing tri-bimaximal mixing for 
neutrinos with the accuracy required by the data, also leads to a natural description 
of the observed pattern of quark masses and mixings. 
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1 Introduction 



It is an experimental fact [Ij that within measurement errors the observed neutrino mixing 
matrix is compatible with the so called tri-bimaximal (TB) form, introduced by Harrison, 
Perkins and Scott (HPS) [2]. The best measured neutrino mixing angle 612 is just about 
la below the HPS value tan^ 612 = 1/2, while the other two angles are well inside the la 
interval [1]. In a series of papers [SJ HI [5l EJ [7j it has been pointed out that a broken flavour 
symmetry based on the discrete group A4 appears to be particularly suitable to reproduce 
this specific mixing pattern as a first approximation. Other solutions based on alternative 
discrete or continuous flavour groups have also been considered [8l[9], but the models 
have a very economical and attractive structure, e.g. in terms of group representations 
and of field content. In most of the models A4 is accompanied by additional symmetries, 
either continuous like U(l) or discrete like Z^r, which are necessary to eliminate unwanted 
couplings, to ensure the needed vacuum alignment and to reproduce the observed mass 
hierarchies. In this way one can construct natural models where the corrections to TB 
mixing can be evaluated in a well defined expansion. 

Recently much attention has been devoted to the question whether a model for HPS 
mixing in the neutrino sector can be suitably extended to also successfully describe the 
observed pattern of quark mixings and masses and whether this more complete frame- 
work can be made compatible with (supersymmetric (SUSY)) SU(5) or SO (10) grand 
unification. Early attempts of extending models based on ^4 to quarks [lOl [6j and to 
construct grand unified versions [TT] so far have not been completely satisfactory, e.g. do 
not offer natural mechanisms for mass hierarchies and for the vacuum alignment. A direct 
extension of the model to quarks leads to the identity matrix for Vckm in the lowest 
approximation, which at first looks promising. But the corrections to it turn out to be 
strongly constrained by the leptonic sector, because lepton mixings are nearly TB, and 
are proven to be too small to accommodate the observed quark mixing angles [S]. Also, 
the quark classification adopted in these models is not compatible with ^4 commuting 
with SU(5) 0. Due to this, larger discrete groups are considered for the description of 
quarks and for grand unified versions with approximate TB mixing in the lepton sector. A 
particularly appeahng set of models is based on the discrete group T', the double covering 
group of A4 [13]. In ref. [H] a viable description was obtained, i.e. in the leptonic sector 
the predictions of the A^ model are reproduced, while the T' symmetry plays an essential 
role for reproducing the pattern of quark mixing. But, again, the classification adopted 
in this model is not compatible with grand unification. Unified models based on the dis- 
crete groups T' [15], 5*4 [16] and A(27) [17J have been discussed. Several models using the 
smallest non-abelian symmetry S3 (which is isomorphic to D^) can also be found in the 
recent literature [TB] . 

In conclusion, the group ^4 is considered by most authors to be too limited to also 
describe quarks and to lead to a grand unified description. In the present work we show 
that this negative attitude is not justified and that it is actually possible to construct a 
viable model based on A4 which leads to a grand unified theory (GUT) of quarks and 

^In ref. [12] an A4 model compatible with the Pati-Salam group SU(4)x SU(2)iX SU(2)fl has been 
presented. 



leptons with TB mixing for leptons. At the same time our model offers an example of 
an extra dimensional GUT in which a description of all fermion masses and mixings is 
attempted. The model is natural, since most of the small parameters in the observed 
pattern of masses and mixings as well as the necessary vacuum alignment are justified by 
the symmetries of the model. For this, it is sufficient to enlarge the flavour symmetry 
by adding a U(l) of the Froggatt-Nielsen type and to suitably modify and extend the 
classification under the flavour group so that finally all fermions transform in an SU(5) 
compatible way. In addition, a Z3 symmetry must be assigned to the fields of the model 
which is, however, flavour-independent. The formulation of SU(5) in extra dimensions 
has the usual advantages of avoiding large Higgs representations to break SU(5) and of 
solving the doublet-triplet sphtting problem. A further ingredient of the model is a U(l)/j 
symmetry which contains the discrete i?-parity as a subgroup. A see-saw realization in 
terms of an A4 triplet of right-handed neutrinos N ensures the correct ratio of light neutrino 
masses with respect to the GUT scale. In the present model extra dimensional effects 
directly contribute to determine the flavour pattern, in that the two lightest tenplets Ti 
and T2 are in the bulk (with a doubling Tj and T/, « = 1, 2 to ensure the correct zero mode 
spectrum), whereas the pentaplets F and T3 are on the brane. The hierarchy of quark 
and charged lepton masses and of quark mixings is determined by a combination of extra 
dimensional suppression factors for the first two generations and of the U(l) charges, while 
the neutrino mixing angles derive from A4. The choice of the transformation properties of 
the two Higgses H^, and is also crucial. They are chosen to transform as two different 

singlets 1 and 1'. As a consequence, mass terms for the Higgs colour triplets are not 
directly allowed |§ and their masses are introduced by orbifolding, a la Kawamura [T^ . 
Finally, in this model, proton decay is dominated by gauge vector boson exchange giving 
rise to dimension six operators. Given the relatively large theoretical uncertainties, the 
decay rate is within the present experimental limits. 

The resulting model is shown to be directly compatible with approximate TB mixing 
for leptons as well as with a realistic pattern of fermion masses and of quark mixings in a 
SUSY SU(5) framework. 

2 The Model 

We consider a SUSY GUT based on SU(5) in 4-1-1 dimensions. Leaving aside extra dimen- 
sional effects for a moment, from the four-dimensional (4D) point of view matter fields are 
chiral supermultiplets transforming as 10, 5 and 1 under SU(5). Part of the flavour sym- 
metry is related to the discrete group A4, whose properties are summarized, for instance, 
in section 2 of ref. 0, whose conventions are adopted here. The three 5 and the three 
singlets (corresponding to the right-handed neutrinos) are grouped into A4 triplets F and 
N, while the tenplets Ti, T2 and T3 are assigned to 1", 1' and 1 singlets of A^, respectively 
(see table 1). The Higgs chiral supermultiplets that break the electroweak symmetry are 
and ifg, transforming as (5,1) and (5,1') under SU(5)xA4. We also consider a set 
of flavon supermultiplets, all invariant under SU(5), that break the A4 symmetry: two 

^Even after A4 breaking they are forbidden at all orders by the U(l)ij symmetry. 



triplets ^pt and Lps and two singlets ^ and ^. The alignment of their vacuum expectation 
values (VEVs) along appropriate directions in flavour space will be the source of TB lepton 
mixing. It is well-known that, for this to work, each triplet should mainly contribute to 
the mass generation of a specific sector. At the leading order and after spontaneous ^4 
breaking, ipsi ^ and ^ should give mass to neutrinos only, while lpt gives mass to charged 
leptons and to down quarks. This separation can be realized with the help of an additional 
spontaneously broken Z3 symmetry under which iV, F, Tj, -ffs 5, (ps-, ^ and ^ are multiplied 
by = exp(227r/3), while (fx is left invariant. The generation of the up quark masses as 
well as the quark mixings will be discussed below. 

The breaking of the grand unified symmetry is a potential source of serious problems, 
like those related to the doublet-triplet splitting and to proton decay. One of the most 
efficient mechanisms to break SU(5) and avoid these problems is the one based on com- 
pactification of extra spatial dimensions [19j. The simplest setting is an SU(5) gauge 
invariant five- dimensional (5D) theory where the fifth dimension is compactified on a cir- 
cle of radius R. The gauge fields, living in the whole 5D space-time, are assumed 
to be periodic along the extra dimension only up to a discrete parity transformation VL 
such that the gauge fields of the SU(3) xSU(2) xU(l) subgroup are periodic, while those 
of the coset SU(5)/SU(3)xSU(2)xU(l) are antiperiodic. Only the gauge vector bosons 
of SU(3)xSU(2)xU(l) possess a zero mode. Those of SU(5)/SU(3)xSU(2)xU(l) form 
a Kaluza-Klein tower starting at the mass level 1/R. From the viewpoint of a 4D ob- 
server, these boundary conditions effectively break SU(5) down to the Standard Model 
(SM) gauge group, at a GUT scale of order 1/R. The transformation f2 is an automor- 
phism of the SU(5) algebra, so that the whole construction can be carried out within an 
SU(5) invariant formalism. An important advantage of this mechanism is that it provides 
a simple solution to the doublet-triplet splitting problem. The parity VL is consistently 
extended to the Higgs multiplets i^s and ifg, also assumed to live in the whole 5D space, 
in such a way that the electroweak doublets are periodic, whereas the colour triplets are 
antiperiodic. In this way we have zero modes only for the doublets and the lightest colour 
triplets get masses of order 1/R. Notice that, if the model is supersymmetric as in the 
case under discussion here, we have an effective 4D N = 2 SUSY, induced by the original 
= 1 SUSY in five dimensions. To reduce N = 2 down to A^ = 1 it is convenient to com- 
pactify the fifth dimension on the orbifold / Z2 rather than on the circle S^. The orbifold 
projection eliminates all the zero modes of the extra states belonging to A^ = 2 SUSY and 
also those of the fifth component of the gauge vector bosons. The zero modes we are left 
with are the 4D gauge bosons of the SM, two electroweak doublets and their A^ = 1 SUSY 
partners. To complete the solution of the doublet-triplet splitting problem, we should also 
forbid a large mass term H^H^, which would otherwise lift the doublet masses. As will be 
explained below, this is automatically guaranteed by the U(l)ij symmetry that we specify 
in table 1. 

For the gauge vector bosons and the Higgses and ifg we will adopt this setup, 
which is described in detail in refs. [20]. For the remaining fields we have much more 
freedom [201 EI]- Indeed the orbifold /Z2 corresponds to a segment where the fifth 
coordinate y runs from to vri?. The boundaries of the segment determine two 4D slices 



of the original 5D space-time. When boundary conditions are consistently defined for the 
local parameters of SU(5) gauge transformations, we find that such transformations are 
generically non- vanishing only in the bulk and at y = 0. At the opposite endpoint of the 
segment, y = ttR, the only gauge transformations that are different from zero are those of 
the SM. Therefore we have three qualitatively different possible locations for the remaining 
fields: in the bulk, at the SU(5) preserving brane y = 0, or at the SU(5) breaking brane 
y = 7tR. We choose to put the two tenplets Ti and T2 of the first and second family in the 
bulk. As explained in ref. [20l |2T] to obtain the correct zero mode spectrum with intrinsic 
parities compatible with symmetry and orbifolding, one must introduce two copies of each 
multiplet with opposite parity fl in the bulk. Therefore Ti 2 is a short notation for the 
copies Ti 2 and T{ 2- The zero modes of Ti 2 are the SU(2) quark doublets Qi,2, while those 
of are and £^2- All remaining = 1 supermultiplets are assigned to the SU(5) 
preserving brane at y = 0. 
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Table 1: Fields and their transformation properties under SU(5), A4, U(l), and U(l)/j. 
Ti and T2 come in two replicas with the same quantum numbers, except for the intrinsic 
parity Q. For simplicity, we only show one of them in the table. 



An interesting feature of the 5D setup is the automatic suppression of the Yukawa 
couplings for the fields living in the bulk. Indeed, a bulk field B and its zero mode B^ are 
related by: 

B = ^=B'' + ... (1) 

where dots stand for the higher modes. This expansion produces a suppression factor 

1 



VttRA 



<1 . (2) 



Thereby, A denotes the ultraviolet cut-off. Such a suppression factor enters the Yukawa 
couplings depending on the field B^. As a result, the hierarchies among the charged fermion 
masses are partly due to the geometrical dilution of the Yukawa couplings involving Ti^2- 
However this dilution cannot account for all the observed hierarchies and, to achieve a 
realistic mass spectrum, we also exploit the Froggatt-Nielsen mechanism. The tenplets Ti 
and T2 are charged under a U(l) fiavour group, spontaneously broken by the VEVs of two 
fields 6 and 6" both carrying U(l) charges —1. The elements of the charged fermion mass 
matrices are provided by higher-dimensional operators with powers of 6 and 6" balancing 



the U(l) charge of the relevant combination of matter fields. Indeed, we need two fields, 9 
and 9" , in order to reproduce a realistic pattern of quark masses and mixing angles. Under 
yl4, 9 is invariant, while 9" transforms as 1". All this is summarized in table 1. 

Notice that, once we have introduced all the fields with the quantum numbers dis- 
played in table 1, there will be no contribution coming from colour triplet exchange to 
the dangerous dimension five operator that induces proton decay in SUSY theories. Actu- 
ally that operator is strictly forbidden as long as the U(l)ii; symmetry remains unbroken. 
Indeed, the superpotential of the effective = 1 SUSY should have U(l)ij charge +2, 
to compensate the i?-charge —2 coming from the Grassmann integration measure d?9. 
With the R assignment in table 1, all superpotential couplings bilinear in the matter fields 
A^, F and T have i?-charge +2 and are allowed. At the same time dangerous operators 
are forbidden. First of all these include the mass term H^H^, that would spoil the solu- 
tion to the doublet-triplet splitting problem. Moreover, since U(l)/j contains the discrete 
-R-parity, also all renormalizable baryon and lepton number violating operators, such as 
FH^ and FFT, are not allowed. Finally, the dimension five operator FTTT, leading to 
proton decay, has i?-charge +4 and therefore is absent. As discussed in detail in ref. [6] 
and briefly recalled in section 4, the U(l)/j symmetry plays also an important role in the 
dynamics that selects the correct vacuum of the theory, which is a crucial feature to re- 
produce nearly TB mixing in the lepton sector. The U(l)ij symmetry is a remnant of the 
SU(2)/j symmetry of the N = 2 SUSY bulk action, before compactification. By reducing 
N = 2 down to A^ = 1 the orbifold projection breaks SU(2)/j down to Eventually, 
after the inclusion of A^ = 1 SUSY breaking effects, the U{1)r symmetry will be broken 
down to the discrete i?-parity, at the low energy scale msusY- The operator FTTT might 
be generated, but with a highly suppressed coupling of the kind {msusY / / n > 0. 
Therefore, the leading contribution to proton decay comes from gauge vector boson ex- 
change and the corresponding proton decay rate is typically small enough, though suffering 
from considerable uncertainties 122]. 



3 Fermion Masses 

The N = 2 SUSY invariance is broken down to A^ = 1 by the orbifold projection, but it 
still forbids 5D superpotential couplings. These couplings should be strictly localized at 
one of the two branes. By choosing the brane at y = 0, the brane action reads: 

d'^x J dy J d^9 w{x)5{y) + h.c. = J d'^x J d^9 w{x) + h.c. . (3) 

The superpotential w, which can be expressed in terms of A^ = 1 superfields, can be 
decomposed into several parts: 

W = Wup + Wdown + Wv + Wd + -- ■ (4) 

The first three contributions in eq. (jl]) give rise to fermion masses after A^, U(l) and 
electroweak symmetry breaking. They are of the form: 

1 9" 9"^ 99"^ 

Wup — J^H^T^T^ + —H^T2T^ + -jjj7^1i^T{T2 + -j^E^T{T?^ 



96"'^ 9^9" 9^9"^ 

Wd^n = ^H-,{FipTy% + -^H-,{FipT)% + ^H-,{FipT)T, + ^H-,{^^^ 

nil 9^9" 99"^ 

+ —H-^{F^t)"T2 + ^H-^{F^t)'Ti + ^H-^{F^t)"T, + ... , (6) 

where dots stand for higher-dimensional operators. In both, Wup and Wdown, the dimension- 
less coefficients of each independent operator have been omitted, for notational simplicity. 
They are not predicted by the flavour symmetry, though they are all expected to be of 
the same order. The powers of the cut-off A are determined by the dimensionality of the 
various operators, by recalling that brane and bulk superfields have mass dimensions 1 and 
3/2, respectively. Some combinations of matter fields, as for instance T1T2 in Wup, appear 
several times, but with the same cut-off suppression. Provided 9 and 9" develop VEVs of 
similar size, the corresponding contributions to the charged fermion mass matrices will be 
of the same order. The bulk matter supermultiplets Ti and T2 come in two copies and, 
to keep our notation compact, the previous formulae do not contain all possible terms 
originating from such a doubling. For instance, F1T2 stands for both combinations F1T2 
and -F1T2, which are suppressed by the same power of A, but can differ by order-one rela- 
tive weights. It is important to keep this point in mind, since it allows to escape the too 
rigid mass relations between the first two generations of charged leptons and down quarks 
predicted by the minimal SU(5) GUT. 

Neutrinos have both Dirac and Majorana mass terms, induced by: 

= j^H,{NF) + {xai + iai){NN) + x,{^sNN) , (7) 

where ^ is defined as the combination of the two independent ^-type fields which has a 
vanishing VEV. Therefore, it does not contribute to the neutrino masses. 

The last term in eq. (jl]), Wd-, is responsible for the alignment of the fiavon fields ipxi 
ips, ^ and C,- The fields 9 and 9" get VEVs from the minimisation of the D-term of the 
scalar potential. We will discuss these issues in the next section. For the time being we 
assume that the scalar components of the supermultiplets acquire VEVs according to the 
following scheme: 

%^ = (..,0,0) , ^ = ivs,vs,vs) , f = - , 

The Higgs multiplets live in the bulk and what matters for the Yukawa couplings are the 
values of the VEVs at ?/ = 0: 

{H,{0)) = -^ , {Hm) = ^ , (9) 

VTrit Vvrit 



where v^^ have mass dimension 1. The electroweak scale is determined by the relation: 
.^^(174GeVf , vl^ rdy\{H,{y))f , vl^rdym{y))f .(10) 

JO JO 



Notice that the electroweak gauge boson masses depend on the 5D averages of [(i^s §(?/)) 
rather than on the values ai y = 0. If the VEVs of H^ ^ are constant along the fifth 
dimension, then = Vu and = Vd- However, if the profile of {H^i{y)) is not fiat in 
the parameters v^^^l less constrained. In order to obtain v^^^l '^u,d, we need some 
special dynamics on the y = and y = ttR branes, that we cannot control without detailing 
additional features of the model, such as the breaking of the residual = 1 SUSY and 
the generation of a non-trivial potential for the electroweak doublets. In this section we 
consider u° ,^ ^ Vu,d as an open possibility and we will discuss a possible application of it. 
All the other fields have vanishing VEVs. 

From these VEVs, the superpotential terms in eqs. flHBllTl) and the volume suppression 
s of eq. ([2]), it is immediate to derive the fermion mass matrices. In the up and down 
quark sector we get, up to unknown coefficients of order one for each matrix element and 
by adopting the convention fRmffi'- 
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where the dots stand for subleading contributions, that will be fully discussed in section 5. 
Here we explicitly see the interplay between the volume dilution and the Froggatt-Nielsen 
mechanism, to achieve the hierarchical pattern of the quark mass matrices. Realistic values 
of quark mass ratios and mixing angles are obtained by assuming 
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We anticipate that, in the absence of corrections to the vacuum alignment given in eq. ([8]), 
the dots receive contributions from highly suppressed operators. In this case the entries 
12, 13 and 23 of m^/ (vtv^) would be of order A^, A^ and A^, respectively. Since vt ~ O(A^) 
(see below), rrib/mt ~ vtv^/v^ ~ A^ is easily reproduced by taking f° ^ Notice that 
there is an overall factor s ~ 0(A), coming from the normalization of the Higgs VEVs, 
eq. ([9]), suppressing both m„ and nid. In order to avoid large dimensionless coefficients, 
we make use of the freedom related to the boundary values f ° ^ and we will assume that 
Vu,d ~ Af ° ^. In this way, the Yukawa coupling of the top quark is of order one and, by 



the patterns given in eqs. (I14fl5p . also all the other couplings are of the same order. 
Alternatively, if the Higgs VEVs are flat along the fifth dimension and f ° ^ = Vu,d, we 
must assume that all Yukawa operators in w have similar couplings of order 1/A [23]. To 
correctly reproduce the quark mixing angle between the first and the second generation, 
a moderate tuning is needed in order to enhance the individual contributions from the up 
and down sectors, which are both of order A^. 

The mass matrix for the charged lepton sector is of the type: 
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mj is relaxed. Indeed, while the third 



We observe that the minimal SU(5) relation 
column of exactly coincides with the third row of rrie, thus implying rrib ~ rrir at the 
GUT scale, the remaining entries are only equal (up to a transposition) at the level of the 
orders of magnitude, since Ti 2 are doubled. This allows to evade the too rigid relations 
rrifj, = rris and rrie = rrid of minimal SU(5). In our 5D setup these relations hold only 
up to order one coefficients and acceptable values of the masses for e, /i, d and s can be 
accommodated. 

In the neutrino sector, after the fields ips and ^ develop their VEVs, the gauge singlets 
N become heavy and the see-saw mechanism takes place. The mass matrix for light 
neutrinos is given by: 
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The neutrino mass matrix is diagonalized by the transformation: 



(17) 



U'^rriyU = diag(mi, m2, ms) 



(19) 



where, in units of s (f°) /A, 
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and U is given by 
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Note that, in the leading approximation, the model predicts the relation: 



(20) 
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It is expected to hold up to corrections of O(A^), as will be discussed in section 5. Notice, 
that in our conventions mj (i = 1,2,3) are in general complex numbers, so that the 
previous relation cannot be used to exactly predict one physical neutrino mass in terms 
of the other two ones. Nevertheless, it provides a non-trivial constraint that the neutrino 
masses should obey. 

To get the right solar mixing angle, we should impose \m2\ > \mi\ and this requires 
COS0 > —\z\/2, where z — h/a and is the phase difference between the complex numbers 
a and h. The neutrino spectrum can have either normal or inverted mass ordering. If 
max(— 1, — 12;|/2) < cos0 < the ordering is inverted, \m^\ < |mi| < |m2|, while \z\/2 < 
COS0 < 1 gives rise to a normal ordering, |mi| < \m2\ < l^^sl- By defining 

r = ^T^loil ^i^ltm > Am^„; = |m2|^-|mi|^ , Am^^^ = | Irngj^ - |mi|^| , (23) 
we find 

.4 . (24) 

We see that a sufficiently small r requires z not to far from either +1 (cos0 = 1, normal 
hierarchy) or -2 (cos0 = —1, inverted hierarchy). If we expand z around +1, we obtain: 

1^1 r = ^^"^atm^+- 
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\m2\ = -Am„i„ r 



ImeeP = ^Am^,^r + ... , (25) 

where we have expressed the parameters in terms of Am^^^ and r. Dots denote terms of 
order and |mee| is the effective mass combination controlling the violation of the total 
lepton number in neutrinoless double beta decay. It is useful to estimate the cut-off A. 
We have roughly 



^ml,^^-^ . (26) 

By taking y^Am^ = 0.05 eV, s^{vlf = (100 GcV)^ and ^ ^ 0.2, we obtain |a|A ^ 10^^ 
GeV, not far from the unification scale. For m ^ f 5 pa the cut-off A is then above 10^^ 
GeV. If we expand z around -2, we get: 

1"^^!' = G + T2)^"^«- + - 

, ,2 /9 13 \ ^ 2 

^2! = + — rj Am„i^ + ... 



12 
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We now have 



^f^m^"-J^ ■ (28) 



By repeating the previous estimate, we find |a|A ^ 10^^ GeV and A slightly below 10^^ 
GeV. 

Several remarks should be made: 

Concerning the lepton mixing, this is dominated by U, eq. (I2T1) . The contribution 
from the charged lepton sector depends on the entries denoted by the dots in rrie- Putting 
all the dots to zero, the charged leptons affect the lepton mixing through rotations of 
order A^, and A^ in the 12, 13 and 23 sectors, respectively. Operators of dimensions 
higher than the ones, considered so far, are strongly suppressed and provide contributions 
of order A^ to the mixing matrix. These are negligible, since the leading effect comes from 
the modification of the vacuum structure of eq. ([H]), due to higher order terms in the scalar 
potential. We shall discuss this in sections 4 and 5. Eventually, such terms modify only 
slightly the TB mixing pattern. 

Apart from Wjy contributions to neutrino masses and mixing angles might come from 
higher dimensional operators, as for instance 

«™ . (29) 

However, they are completely negligible compared to those discussed above. If we forced 
this type of operator to be the dominant one, by eliminating the singlets N from our 
model, we would need a value of A too small compared with the GUT scale. 

Depending on the value of z, our model gives rise to two separate branches in the 
neutrino spectrum, both characterized by a nearly TB mixing. On the first branch, z ^ +1, 
we find a spectrum with normal hierarchy, while on the second branch, z ^ —2, we get 
an inverted hierarchy. A degenerate spectrum is actually disfavored in our construction, 
since it would require 2; -C 1 (see eq. (1201) ) which leads to r close to 1/2, as can be read 
off from eq. fl2^ . This can obviously not be reconciled with the data. 

In our model the possibility of normal hierarchy is somewhat more natural than the 
one of inverted hierarchy. There is no reason a priori why z should be close to +1 or to 
—2 and reproducing r requires some amount of tuning. However, such a tuning is stronger 
for inverted hierarchy (ih) than for the normal one (nh), as can be seen by 



dr 
dz 



dz 
dr 



. (30) 



nh 

The derivatives are computed at the relevant value of z in each case and r is the exper- 
imental value. Moreover the solution with a normal hierarchy has a domain of validity 
in energy larger by a factor of 1/v^ ~ 5.6 and extends beyond 10^^ GeV. In the normal 
hierarchy solution we find with the help of eq. (!25l) 

^ |mi| ^ (0.06 - 0.07) eV and |mee| ~ 0.007 eV . (31) 

i 

It is interesting to see that \mee\ is close to the upper limit of the range expected in 
the normal hierarchy case, being not too far from the aimed for sensitivity of the next 
generation of neutrinoless double beta decay experiments, 0.01 eV. This is partly attributed 
to the fact that \mi \ 0.005 is different from zero and in part to the absence of a negative 
interference with the 777,3 contribution, as = 0. 



4 Vacuum Alignment 



Here we discuss the minimisation of the scalar potential, in order to justify the VEVs 
assumed in the previous section. We work in the limit of exact SUSY. This will not allow 
us to analyse the electroweak symmetry breaking induced by and ifg, whose VEVs are 
assumed to vanish in first approximation. Indeed all the VEVs we are interested in here, 
i.e. those of the flavon fields ^s,T) ^, C,, and 9", are relatively close in magnitude to the 
cut-off A and therefore much larger than the electroweak scale, which will be consistently 
neglected. Moreover we work at leading order in the parameter 1/A, that is we keep 
only the lowest dimensional operators in the superpotential shown in the previous section. 
Subleading effects will be discussed later on. All the multiplets but the flavon ones are 
assumed to have vanishing VEVs and set to zero for the present discussion. We regard 
the U(l) Froggatt-Nielsen flavour symmetry as local. Since the field content displayed in 
table 1 is anomalous under the U(l), we need additional chiral multiplets to cancel the 
anomaly. These multiplets can be chosen vector-like with respect to SU(5), so that they 
only contribute to the U(l) anomaly. Here we do not need to specify these fields, but we 
must presume that they do not acquire a VEV. Within these assumptions the relevant 
part of the scalar potential of the model is given by the sum of the F-terms and of a 
D-term: 

V^VF + Vn , (32) 



where ipi stands for the generic chiral multiplet. Only the last term in eq. (jl]), Wd, 
contributes to the VEVs we are looking for. It is given by: 

Since also the terms in Wd have to have i?-charge -|-2, we introduce additional gauge 
singlets, so called driving fields, ^Jq, yjf and with i?-charge +2 (see table 1). Note 
that therefore all terms in Wd are linear in these fields. Note further that due to U(l) 
invariance neither the multiplet 6', nor the multiplet 9" is contained in Wd- Moreover the 
D-term Vd does not depend on v^s.t, ^5 ^5 which are all singlets under the (gauged) U(l). 
The expression of Wd and the minimisation procedure are exactly as described in ref. [5] 
and leads to the result anticipated in the previous section: 

((^t) = {vt. 0, 0)A , wtA 



2g 
'4 

— -i 
35-3 



(0 = «A , 

(0 = (34) 

with u undetermined and = 3g^ , (74 = — In the following we take vt, vs and u to 
be of O(A^). This order of magnitude is indicated by the observed ratio of up and down 



or charged lepton masses, by the scale of the hght neutrino masses and is also compatible 
with the bounds on the deviations from TB mixing for leptons. 
The D-term is given byjfl 

Vd = \{Ml, - gFM\e\^ - gFN\e"\^ + ...)' (35) 

where qfn is the gauge coupling constant of U(l) and Mp^ denotes the contribution of 
the Fayet-Iliopoulos term. We have omitted the SU(5) contribution to the D-term, whose 
VEV is zero. There are SUSY minima such that Vp = Vo = 0. The vanishing of Vo 
requires 

QenW + gFN\e"\^ = Ml, . (36) 

If the parameter Mpj is positive, the above condition determines a non-vanishing VEV for 
a combination of 9 and 9" . Here we assume that the VEVs fulfil t, t" ~ 0(A) according 
to eqs. (PITSjl . The different order of t, t" versus vt, vs and u can be attributed to the 
different couplings and mass parameters in Vd and Vp- 

Finally, we discuss the subleading corrections to the vacuum alignment. As already 
noticed above, the fields 9 and 9" cannot couple to the flavon fields, since the flavons ifT-, 
¥^3, ^5 V^o ' V'o ^0 are not charged under the U(l) symmetry, responsible for the 
charged fermion mass hierarchy. Therefore, the subleading effects in the potential arise 
from terms made up of one driving field and three fields (fp, ^ and ^. They induce 
shifts in the VEVs shown above and thereby influence the mass matrices, as discussed in 
the next section. Since the fiavon field content of this model is essentially the same as the 
one in ref. [6], not only the renormalizable part of Wd coincides, but also the subleading 
terms are the same. Hence, we do not need to repeat this discussion and we only state 
the results found there. The shifted VEVs are 

(V3r)/A = {vt + 5vti,5vt2,5vt^) , 

(V95)/A = {vs + 6vi,vs + 6v2,vs + 5v:i) , 

(0/A = u , 

(0/A = 5u' , (37) 

where u remains undetermined and, once we have taken wt,5, u ~ O(A^), all shifts are 
suppressed by a factor of order A^: 5v/v ~ O(A^). As found in ref. [6| the following 
relation holds: 

5vt2 = ^vt3 ■ (38) 

Higher order corrections to t and t" simply amount to a rescaling that does not change 
their individual order of magnitude which remains of 0(A). 



5 Subleading Corrections 

In this section, we analyse the effects of the subleading corrections in terms of A to the 
fermion masses and mixings. The corrections arise from additional insertions of the fiavons 
ipT, 'Ps, ^ and ^ as well as from shifts of the VEVs shown above. 

^Note that |6'"p is a singlet under A4, because 9" ~ 1" and 6"* ^ 1' under A4. 



5.1 Corrections to Wup 

In the up quark sector the leading order terms only involve the fields 6 and 6", since they 
are the only fields which have a non-vanishing U(l) charge among the gauge singlets of 
the model. The subleading terms then additionally involve the fields ipT, V^s, ^ and ^. 
As the tenplets transform as singlets under A4 and the combinations TiTjH^O^'O"^ are 
invariant under the Z3 group, we cannot multiply the Wup terms by a single fiavon field. 
The most economic possibility is to insert two fiavons, namely (pr ^t- Among the three 
contractions leading to a 1 or 1' or 1" representation of A/^ only the 1 has a non-vanishing 
VEV, given that (v^t) = {vt, 0, 0) A. Therefore the dominant subleading corrections to 
the up quark mass matrix have the same structure as the leading order results and are 
suppressed by an overall factor t>|, ~ O(A^). The fields ips and ^, ^ can only couple at 
the level of three fiavon insertions due to the requirement of invariance. However, all 
contributions stemming from three fiavon insertions are suppressed by relative to the 
leading order term. Similarly, the corrections due to shifts in the VEVs contribute at most 
at relative order A^. For the up quark masses and the mixings all these corrections are 
negligible. 



5.2 Corrections to tc^, 



own 



In the down sector the main effect of the subleading corrections is to fill the zeros indicated 
by dots in the upper triangle of m^. In order to maintain the A/^ invariance the leading 
order terms include one insertion of the flavon LfT- The subleading corrections arise from 
two effects: a.) replacing ipx with products of fiavon fields and h.) including the corrections 
to the VEVs of ipx- The replacement of (pT with a product PtPt is the simplest choice 
compatible with the Z3 charges. Note that this is similar to the up quark sector. If the 
VEVs are unchanged this contribution to is of the same form as displayed in eq. fll2p 
and suppressed by vt ~ O(A^) compared to the leading result due to the additional fiavon 
field. Therefore this type of correction does not fill the zeros in m^. They are filled by the 
corrections coming from the VEV shifts inserted in the terms containing one fiavon LfT- 
Considering that we assumed all 5v/v ~ O(A^), the corrections to the matrix elements of 
rud are of the following order in A: 



5mA 



( A^ A^ A2 \ 

A^ A^ A2 
\ A^ A^ A^ j 



vt\v^^ . (39) 



As said, the matrix elements which are already non-vanishing at the leading order, eq. 
(fT2l) . receive additional corrections from the two fiavon insertion ^pt'-Pt- These are of the 
same order as the corrections from the VEV shifts, e.g. for the element 11 also of order 
A^. In summary, the zeroes in the elements 12, 13 and 23 of m^, appearing at leading 
order, are replaced by terms of order A^, A^ and A^, respectively, in units of vt^v^- 

In our model the relation nid = is not valid for the first two families but it still 
holds at the level of orders of magnitude for each entry. So the powers of A are also the 
same for each matrix element of m^m^ and of m^mg. This is important as the matrix 



mime is diagonalized by the unitary matrix Ue that enters in determining the leptonic 
mixing matrix U = U^U^- The results just described for the subleading corrections on m^ 
and ml imply that Ue induces corrections of O(A^) on all mixing angles in U, that is, in 
our case, corrections of O(A^) to the TB values of each mixing angle. 

5.3 Corrections to w,^ 

Also the Wi, term of the superpotential, eq. (JTj), is modified by terms with more flavon 
factors and by subleading corrections to the VEVs. The Dirac mass term, proportional 
to H^{NF), is mainly modified by a single ipT insertion, that produces corrective terms 
suppressed by a O(A^) factor. These corrections are of the same order as those arising for 
Majorana mass terms. In fact, NN can be in a 1 , 1', 1" or 3^ combination. Since A^A^ ~ u'^ 
under Z^, the singlet 1 can be multiplied by ^ (the singlet leading term) or by {(fTfs) 
(which can be absorbed into a redefinition of the leading term), 1' by ((pTfs)'\ 1" by 
{fT^sY and 3s by (ps (the triplet leading term) or by {(frO or {ipT(ps)3s {fT^s)3a- 
two flavon insertions lead to corrections of relative order of O(A^) to the matrix elements 
of the Majorana matrix. In addition, the shifts of the ips VEVs applied to the triplet 
leading term also produce O(A^) corrective terms. As it is easy to check, in general there 
are enough parameters so that all 6 independent entries of the (symmetric) Majorana mass 
matrix receive a different correction at O(A^). 

The described corrections affect the neutrino masses and, together with the corrections 
to me, also all lepton mixing angles. To be compatible with the data, given the accuracy 
of the TB approximation, the dominant corrections must be of O(A^) at most, and this is 
precisely the magnitude of the terms that we have just mentioned. 

6 Conclusion 

We have constructed a SUSY SU(5) grand unified model which includes the description 
of TB mixing for leptons. For this it is not only necessary to adopt an classification 
of quarks and leptons compatible with SU(5), but also to introduce additional U(l) and 
Zn symmetries and to suitably formulate the grand unification model. We find that the 
most attractive solution to cope with the different requirements from fermion mass and 
mixing hierarchies, from the problem of doublet-triplet splitting in the Higgs sector, from 
proton decay bounds and from maintaining bottom tau unification only, is a formulation 
in 5 space-time dimensions with a particular location of the different fields, with some of 
them on the brane at |/ = and some in the bulk. The latter include the gauge and Higgs 
fields as well as the tenplets of the first two, i.e. lightest, families. The resulting model 
naturally leads to TB mixing in first approximation with corrections of O(A^) from higher 
dimensional effective operators, together with reproducing the observed mass hierarchies 
for quarks and charged leptons and the CKM mixing pattern. In the quark sector, however, 
as is typical of U(l) models, only orders of magnitude are determined in terms of powers of 
A with exponents fixed by the charges. A moderate fine tuning is only needed to enhance 
the CKM mixing angle between the first two generations, which would generically be of 



0(A2), and to suppress the value of r, given in eq. (12^ . which would typically be of order 
1. The latter feature is also true in all purely leptonic A4 models, in which ^4 leads to 
the correct mixing, but not directly to the spectrum of the neutrino masses. Actually the 
model allows for both types of neutrino mass hierarchy, the normal and the inverted one. 
The normal hierarchy is, however, somewhat more natural, since it requires less tuning 
to reproduce r. Furthermore, it is consistent with a larger value of the cut-off A. If the 
normal hierarchy is the correct one, the model predicts the sum of neutrino masses to be 
around (0.06 — 0.07) eV and |mee| to be close to 0.007 eV. Therefore, \mee\ is not far from 
the future aimed for experimental sensitivity. Finally, all subleading corrections to the 
leading order result of fermion masses and mixings have been carefully analysed. 

In conclusion, we have demonstrated that the simple A4 approach to TB mixing is 
compatible with a grand unified picture describing all quark and lepton masses and mix- 
ings. 
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